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The mode-converted ion-Bernstein wave excited in tokamaks is shown to exhibit certain very 
interesting behavior, including the attainment of very small poloidal phase velocities, the reversal 
of poloidal direction, and up-down asymmetries in propagation and damping. Because of these 
effects, this wave holds promise for channeling a-particle power to ions, something that would make 
a tokamak fusion reactor far more attractive than presently envisioned. 
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One way to achieve significantly improved tokamak 
performance is to direct fusion a-particle power into ei- 
ther electrons to sustain the plasma current [1] or into 
ions to increase the plasma reactivity [2]. Such a chan- 
neling of the a-particle power could make possible the 
operation of fusion reactors in regimes in which there is 
a significant hot, nonmaxwellian component to the fuel 
ions, or in which the fuel ion temperature can be much 
greater than the electron temperature [3,4]. The fusion 
reactivity is then greatly enhanced, making possible a 
much more economical reactor. 

The free energy in the a-particles, namely, the free en- 
ergy associated with the gradient in a-particle pressure, 
may be tapped through a resonant quasilinear interaction 
with a poloidally propagating short wavelength electro- 
static wave [1]. Such a wave diffuses the a-particles both 
in space and energy, rather than just in energy. The free 
energy may be accessed if k±p a ^> 1, where p a = v a /Q, a 
is the a-particle gyroradius, and where, for radial dif- 
fusion, ki_ ~ kg. Upon interaction with the wave, the 
particle velocity is assumed displaced instantaneously in 
the 6 direciton by Ai>j_, in the paralle direction by Any, 
and as a result the position of the guiding center changes 
by Av A i||/fi a . (The confining magnetic field defines 
the || and _L directions; in cylindrical coordinates, 6 is 
the azimuthal or poloidal direction.) 

Importantly, the change in the radial guiding center 
displacement, Ar gc = Avj_/£l a , is related directly to the 
energy change, AE = m a (vgAv±_ + i>||Ai>||). This occurs 
because particles in a magnetic field are diffused by waves 
along the paths 



v±Av±_ + — ui/k\\)Av\\ = 0, 



(1) 



and, since k±p a ^> 1, an a-particle is unmagnetized in 
the wave, and, hence, Landau resonant if at any point in 
its orbit 



vg = (ui - k^)/kg. 
It follows from Eqs. (1) - (2) that 

Ar gc = AEkg/m a Q a uj. 
The usual resonance condition 

v\\ = (uj - nQ. a )/k\\. 



(2) 



(3) 



(4) 



is recovered by requiring successive kicks add over many 
gyroperiods. 

It is Eq. (3) that leads to the quest for large-fcg waves. 
This is a necessary, if not sufficient, requirement, for a- 
channeling, and the focus of this paper. Clearly, in or- 
der to extract successfully a-particle power, the key ratio 
Ar gc /AE must neither be too large nor too small. In 
the former case, the a-particle leaves the tokamak be- 
fore giving up its energy. In the latter case, it is difficult 



to establish a population inversion. The ideal condition 
is for this ratio to be such that Ar gc is on the order of 
the minor radius for extracting the full a-particle energy. 
From Eq. (4), this maximum excursion is 



Ar n 
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where p a o ~ 4cm is the initial gyroradius for 3.5 MeV 
a-particles. For 1 < Q. a /uj < 3/2, corresponding to 
wave frequencies between the deuterium and tritium res- 
onances (such as the IBW), then for large kg sufficiently 
large, say kg ~ 4 cm -1 , one obtains 32 cm < Ar max < 
48 cm, which is a range suitable for extraction, although 
toroidal effects can be expected to modify somewhat this 
range. The required diversion of power occurs if the wave 
then damps on fuel ions. 

To achieve large kg , one might consider the lower hy- 
brid wave [1], but the wave may be difficult to excite 
in the proper regime, or there may be other difficulties 
[5], The mode-converted ion Bernstein wave [6], however, 
might be excited more easily with the requisite poloidal 
phase velocity and high wave number. This can occur as 
follows [7]: In a D-T (or other multispecies) plasma, the 
ion Bernstein wave (IBW) may be excited by launching 
a fast wave that mode-converts at the ion-ion resonance 
surface into an IBW [8]. This mode-converted ion Bern- 
stein wave propagates essentially in the poloidal direction 
if two choices are made: one, the ion-ion hybrid resonance 
layer, which is more or less a vertical plane, is chosen 
to intersect near the tokamak magnetic axis; and, two, 
the fast wave is launched so as to intersect this vertical 
plane somewhat off the magnetic axis. Thus, horizontal 
is equivalent to poloidal propagation. 

Of course, coupling to an IBW is possible through a 
variety of other methods [9]. One possibility would be to 
launch a slow electrostatic wave (electron plasma wave), 
which will subsequently convert to an IBW at the lower 
hybrid resonance. It would be difficult, however, in such 
a manner, to achieve the high poloidal wave numbers that 
are necessary, and that are achieved here, through mode 
conversion at a vertical ion-ion hybrid resonance. The 
method of coupling analyzed here, in contrast, exploits 
the production of the short wavelengths characteristic of 
a resonance layer that is aligned precisely perpendicular 
to the poloidal direction. 

The location of this ion-ion hybrid resonance is deter- 
mined by [10] 



f D n T + f T n D 
f D n D + f T n T ' 



(6) 



where u> is the wave frequency, Sl^ = eB/mo is the deu- 
terium gyrofrequency, Sly is the tritium gyrofrequency, 
and fo and fx are their respective fractional densities. 
The horizontal location of this resonance, which may be 
chosen so as to intersect the magnetic axis, may be con- 
trolled through the toroidal magnetic field, B, through 
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the wave frequency, u>, or through the relative concentra- 
tion of the different ionic species. 

Note that, since the gradient of the magnetic field 
strength in the major radius direction (R) is the principal 
inhomogeneity intrinsic to the mode conversion process, 
the R component of the wavevector of the outgoing wave 
will increase most rapidly as the IBW propagates away 
from the mode conversion layer, thus achieving essen- 
tially poloidal propagation with high poloidal wavenum- 
ber. 

In our analysis of the propagation of such a wave in 
a reactor-grade plasma, several things have been discov- 
ered: One, through a suitable choice of wave and toka- 
mak parameters, and under conditions [4] important for 
channeling a-particle power, namely ion temperatures 
about 20 KeV, and electron temperatures about 10 KeV, 
high kg can in fact be reached. Two, in this regime, the 
IBW can be made to reverse direction, so that the fields 
near high kg will be enhanced, as will the channeling ef- 
fect. Three, it is possible to augment and exploit certain 
up-down asymmetries (because of the poloidal magnetic 
field) both for off-axis current-drive and for channeling a- 
particle power. Four, a deuterium rich mixture achieves 
simultaneously both the high-kg channeling effect and the 
diverting of power to ions. 

The dispersion relation, which describes the ion Bern- 
stein wave, may be written in the form (see Eq. (1) of 
Ref. [9]) 

where the Kjj are elements of the susceptibility tensor, 
and n = kc/w is the index of refraction. In the regime 
important for channeling, it is not possible [11] to em- 
ploy conventional (small k±p) approximations [12-14] of 
these susceptibilities. In what follows, the wave propa- 
gates horizontally (x) and toroidally (z), but not verti- 
cally (k y = 0). 

Because our interest is primarily in off-axis effects, the 
poloidal field importantly relates the parallel to the per- 
pendicular wavenumber through fey ~ k z + 9k x , where 
9 = Bg/B. The first term, k z , is the launched paral- 
lel wavenumber or the parallel offset. The second term 
represents the projection of the horizontal wavenumber, 
k x ~ kg, on the magnetic field in the presence of a small 
poloidal magnetic field, Bg, which is assumed horizon- 
tal in the regime of interest. Suppose that the antenna 
launching the fast wave is up-down symmetric, and sit- 
uated on the low-field periphery of the tokamak. Then 
one might envision two rays of fast waves propagating 
from the antenna, one that undergoes mode conversion 
to the IBW in the top half of the tokamak and the other 
in the bottom half. The resulting IBW spectrum, how- 
ever, is not symmetric for nonzero poloidal field, since Bg 
is up-down antisymmetric. 



In Fig. 1, we trace such rays after the mode conversion 
for a 50:50 DT mixture. Note that for the case of can- 
cellation of fey (6 < 0, or, say, for the IBW propagating 
in the top half of the tokamak) , there is a tendency for 
the wave to reach sufficiently high k± to be of interest for 
channeling. Note also that, particularly for k z = 0.1, the 
wave in fact turns around prior to damping, returning to 
the mode conversion layer at an upshifted wavenumber 
[15]. A caveat to the results here, and to be explored fur- 
ther, is that the ray tracing approach may not be valid, 
particularly for 9 > 0, where the damping is so intense as 
to occur in a distance short compared to a wavelength. 
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FIG. 1. Trajectory of IBW k x (x) from dispersion rela- 
tion of Eq. (7), with parameters: lo = 1.69 • 10 8 rad sec -1 , 
n e = 10 14 cm -3 , T e = lOkeV, T D = T T = 20keV, with frac- 
tional densities, fo = fr = 0.5. The plasma inhomogeneity in 
the x direction is due entirely to the variation in magnetic field 
B = 4.8 • 10 4 (z + 6>x) R/x gauss, with 6 = ±0.05, and " major 
radius" R = 260 cm. The broken curves are trajectories for 
the values .05, .75, and .1 cm -1 of the toroidal wavevector k z . 
The solid lines connect the coordinates at which .75, .5, and 
.2 of the initial power remains. 

For channeling alpha power, it is significant both that 
the horizontal group velocity passes through zero and 
that this can be made to occur at high kg . There is field 
enhancement where the group velocity vanishes, so that 
less power is required for the waves to diffuse a-particles 
in a time short compared to their slowing down time. 

Note also that there are considerable up-down asym- 
metries. First, clearly, the sense of poloidal propagation 
is opposite. Second, for 9 > 0, high-fc e are not reached, 
since this wave damps immediately. Third, the damping 
for 9 > is on ions, whereas, for 9 < 0, it is on electrons. 
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These asymmetries lead to different but not necessarily 
canceling effects. If the fast wave antenna is poloidally 
phased, then it would be possible to launch just the wave 
that produces the more desirable effect. However, this 
phasing is not easily achieved, and, in any event, present 
tokamaks are generally not equipped for such phasing. 
Hence, it is of interest, therefore, to examine whether, 
for example, extraction of a-particle power can occur by 
the mode-converted IBW with toroidal phasing only. 

For the extraction process, it is clear that the large 
disparity in kg leads to significant movement in the a- 
particles from the top interaction with insignificant can- 
cellation from the bottom one. In addition, there would 
be amplification of the wave by a-particles at the top, 
but damping at the bottom. The only loss in not having 
poloidal phasing is that half the power is wasted with 
respect to extraction. 

The current drive effect in the midplane [16], i.e. for 
9 = 0, should not be affected by poloidal fields. The cur- 
rent drive effect off axis, however, is more complicated. 
Since the wave travels in opposite directions on top and 
bottom, and in both cases is completely damped, the 
current drive effect might be thought to vanish off axis, 
and, in fact, has been shown to vanish at least for very 
small fey [17]. However, for realistic fey, there are inter- 
esting asymmetries: First, in the top case, the wave is 
90% damped on electrons, whereas in the bottom case 
it is 2/3 damped on deuterium, 1/3 on tritium. Hence, 
the electrons and ions absorb momentum in opposite di- 
rections, but their currents from this effect add. Second, 
the parallel velocity of electrons and tritium ions that 
absorb wave power must be of the same sign as w/fcy, 
whereas the parallel velocity of the resonant deuterium 
ions must be of opposite sign. Hence, to the extent that 
the resonant ions absorb power, there is the minority 
species heating current drive effect (MSHCD) that could 
be in opposite directions for these ions [18]. Depending 
on the effective bulk ion charge state, the trapped parti- 
cle fraction, and other parameters, the MSHCD currents 
can actually be made to flow in either direction [19,20]. 
Thus, it is clear that sufficient asymmetry exists for the 
current drive effect even off axis, and that here the effects 
of the oppositely traveling waves can even be additive. 

Note, however, that while extraction of a-particle 
power appears to be in place, and while this power may 
in fact be channeled to produce current drive even with 
toroidal phasing only, it follows from Fig. 1 that this 
power can not be channeled to ions for increased reactiv- 
ity. To accomplish this channeling effect, it will be nec- 
essary to explore further the parameter space. In Fig. 
2, we vary the fuel mix, showing the damping on each 
species for (a) bottom case (i.e. 6 > 0), (b) 9 = 0, and 
(c) top case. (Note that the on-axis case has only 50% 
electron damping, which limits the mode-converted cur- 
rent drive, but the electrons for this case are relatively 
cold.) 
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FIG. 2. Left axis displays fractional power absorbed by 
electrons (long dashes), deuterium (dots), and tritium (short 
dashes). Right axis displays k x for which 50% of power re- 
mains (solid line). All are plotted as a function of fractional 
deuterium density fo, with 9 = .05, 0, —.05 in (a), (b), and 
(c), respectively. 
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FIG. 3. Ray trajectories for fo = -7. Other parameters are 
as in Fig. 1. 

Of particular interest is that, for a 70:30 DT mix, there 
is the opportunity for both high kg and damping on tri- 
tium ions. Here, most of the power can be channeled to 
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tritium ions at near maximum kg . Note that this is quite 
relevant to reactor optimization, since at 70:30, there is 
85% of the maxwellian fusion power, plus increased op- 
portunity for an enhanced reactivity of a nonmaxwellian 
tritium tail [4]. 

In Fig. 3, we show the ray tracing for the 70:30 case, 
which indicates that here, as in the 50:50 case, even 
with toroidal phasing only, the channeling effect remains. 
However, the current drive effect on top and bottom 
largely cancels, since in both cases primarily tritium ab- 
sorbs the power. As in the 50:50 case, the horizontal 
group velocity vanishes, and the dissipation, especially 
in the bottom trajectory, is over a very short distance. 
Also, note that the mode has shifted very much nearer 
to the tritium resonance. 

Since the high-kg now occurs at only u> ~ 0.74 it 
is of interest to check whether there is sufficient Doppler 
shift for a-particles to interact with the wave. From Eq. 
(3), using p a ~ 4cm, and, from Fig. 3, taking k± ~ 
3 cm -1 and fey ~ —.075 cm -1 , we find that v\\/v±_ ~ 1, 
indicating that the resonance still falls within the bulk of 
the a-particle distribution. 

Although the effect reported is neither fully explored 
nor optimized here, it can be concluded from this study 
that large kg waves can be excited in tokamaks, off-axis, 
so as to divert a-particle power to ions. At some sac- 
rifice in optimization, both the diversion of power and 
current drive can be expected, even without poloidally 
asymmetric antenna phasing. 

The channeling effect itself requires testing on toka- 
maks capable of D-T fusion, of which at present there 
are two, JET and TFTR. It would be important, how- 
ever, first to observe the effects predicted here with re- 
spect to the up-down asymmetries in the magnitude of 
the wave vector, energetic particle loss, and current-drive. 
Mainly what is needed to observe these effects is off-axis 
excitation of the ion-Bernstion wave, with added benefit 
if this could be done asymmetrically either poloidally or 
toroidally. 

The authors are indebted to C. F. F. Karney and 
R. Majeski for very useful comments. This work was 
supported by the United States Department of Energy 
under contract number DE-AC02-76-CHO3073. 
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